function [VecResults,VecOutliers2]=RemoveOutliers2(VecTemp,ValueTemp4,ValueTemp5)

VecOutliers2=[];
ValueMean=mean(VecTemp(:,ValueTemp5));
ValueSQRT=sqrt(var(VecTemp(:,ValueTemp5)));
ValueTemp=find(VecTemp(:,ValueTemp5)<ValueMean+ValueTemp4*ValueSQRT&VecTemp(:,ValueTemp5)>ValueMean-ValueTemp4*ValueSQRT);
VecOutliersIndexTemp=ValueTemp;
ValueTemp3=find(VecTemp(:,ValueTemp5)>=ValueMean+ValueTemp4*ValueSQRT);
VecOutliers2=[VecOutliers2;VecTemp(ValueTemp3,ValueTemp5)];
ValueTemp3=find(VecTemp(:,ValueTemp5)<=ValueMean-ValueTemp4*ValueSQRT|VecTemp(:,ValueTemp5)>=ValueMean+ValueTemp4*ValueSQRT);
VecOutliersIndex=ValueTemp3;
VecTemp2=VecTemp;
VecTemp=VecTemp(ValueTemp,:);
while size(ValueTemp3,1)>=1
    ValueMean=nanmean(VecTemp(:,ValueTemp5));
    ValueSQRT=sqrt(nanvar(VecTemp(:,ValueTemp5)));
    ValueTemp=find(VecTemp(:,ValueTemp5)<ValueMean+ValueTemp4*ValueSQRT&VecTemp(:,ValueTemp5)>ValueMean-ValueTemp4*ValueSQRT);
    ValueTemp3=find(VecTemp(:,ValueTemp5)>=ValueMean+ValueTemp4*ValueSQRT);
    VecOutliers2=[VecOutliers2;VecTemp(ValueTemp3,ValueTemp5)];
    ValueTemp3=find(VecTemp(:,ValueTemp5)<=ValueMean-ValueTemp4*ValueSQRT|VecTemp(:,ValueTemp5)>=ValueMean+ValueTemp4*ValueSQRT);
%     VecOutliersIndex=[VecOutliersIndex;VecOutliersIndexTemp(ValueTemp3,1)];
%     VecOutliersIndexTemp=VecOutliersIndexTemp(ValueTemp,:);
    VecTemp=VecTemp(ValueTemp,:);
end
ValueTemp=find(VecTemp2(:,ValueTemp5)<min(VecTemp(:,ValueTemp5)));
VecTemp=[VecTemp;VecTemp2(ValueTemp,:)];
VecResults=VecTemp;